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We investigate effects of dipole-dipole (DD) interactions on immiscibility-miscibility transitions 
(IMTs) in two-component Bose-Einstein condensates (BECs) trapped in the harmonic-oscillator 
(HO) potential, with the components linearly coupled by a resonant electromagnetic field (accord- 
ingly, the components represent two different spin states of the same atom) . The problem is studied 
by means of direct numerical simulations. Different mutual orientations of the dipolar moments in 
the two components are considered. It is shown that, in the binary BEC formed by dipoles with 
the same orientation and equal magnitudes, the IMT cannot be induced by the DD interaction 
alone, being possible only in the presence of the linear coupling between the components, while the 
miscibility threshold is affected by the DD interactions. However, in the binary condensate with the 
two dipolar components polarized in opposite directions, the IMT can be induced without any linear 
coupling. Further, we demonstrate that those miscible and immiscible localized states, formed in 
. the presence of the DD interactions, which are unstable evolve into robust breathers, which tend to 

' keep the original miscibility or immiscibility, respectively. An exception is the case of a very strong 

DD attraction, when narrow stationary modes are destroyed by the instability. The binary BEC 
composed of co-polarized dipoles with different magnitudes are briefly considered too. 
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I. INTRODUCTION 



The development of the trapping techniques has made it possible to create multicomponcnt Bose-Einstein con- 
densates (BEC) formed by atoms in different electronic states The multicomponent BEC is an ideal system for 
studying phase transitions and the coexistence of different phases, which are topics of great significance to many areas 
of physics. The advantage of the gaseous BEC systems is the possibility to very accurately describe their dynamics in 
the framework of the mean-field approximation. This approach does not apply to multicomponent quantum systems in 
. condensed states, because their densities are much higher, making the delta-functional pseudopotential, which is used 
t-H ' to describe the local interactions between atoms in dilute BEC by means of the single scattering-length parameter, 
, inappropriate. 

In particular, experiments have been performed in a binary BEC created in 87 Rb, which contains atoms in two 
different hyperfine (spin) states [l[ . The physics of binary BECs gives rise to novel paradigms for ground-states wave 
function and excitations on top of them [2[ , • One generic scenario concerning the binary condensates leads to an 
equilibrium state characterized by separation of two immiscible species in different spatial domains 0] ■ In addition, a 
• *h , resonant electromagnetic spin-flipping field, with a frequency in the range of several GHz, can induce a linear coupling 
(interconversion) between the two species [5j. This technique opens a way to observe various effects, such as Josephson 
oscillations between the two states [6] , domain walls Q , "co-breathing" oscillation modes Q , nontopological vortices 
@, and others. 

In the absence of the linear coupling between the components, the superfluid components formed by atoms in 
two hyperfine states of 87 Rb, \F = l,mp = —I) and \F = 2, nip = +1), are immiscible, although being close to the 
miscibility threshold. On the contrary, the mixture of \F = l,mp = 1) and \F — 2, rap — 2) states is miscible, but 
the separation between them can be induced by a gradient of the magnetic field [lOj, or by causing them to flow in 
opposite directions . The dynamics of the separation between the two respective superfluid components has been 
studied in detail experimentally [l2j . 

The immiscibility-miscibility transition (IMT) in the binary BEC may be induced through the Feshbach resonance, 
which makes it possible to change the scattering len gth s for inter-atomic collisions by means of an external magnetic 
field [13], or a properly tuned coherent optical signal [14j. In particular, the change of the miscibility of 85 Rb and 87 Rb 
condensates by means of the Feshbach resonance, which affects the intrinsic scattering length of the former species, 
was demonstrated in the experiments [l5j | . 

On the other hand, it has also been shown that the IMT between the components which represent different spin 
states of the same atom may be induced by the linear coupling between them [16[ . In this way, an immiscible binary 
condensate may be converted into a miscible one, using the spin-flipping electromagnetic wave of a rather small 
amplitude. 
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In all the above-mentioned cases, the studies of the binary condensates were limited to the case of contact inter- 
particle interactions. On the other hand, a great deal of attention has been recently drawn to dipolar BECs. They can 
be created in gases of ultracold polar molecules [13] , in chromium [l8| , or in the case when atomic electric dipoles are 
induced by the optical pumping of atoms into a Rydberg state. Binary dipolar BECs are possible too - in particular, 
those formed by atoms in two different Rydberg states [l9j . Another realization of a dipolar mixture is possible in 
condensates of heteronuclear diatomic molecules classified as type A according to the Hund's nomenclature (2^, in 
which the direction of the magnetic moment is correlated (parallel or antiparallel) with the orientation of the molecular 
axis 12111 . Through this mechanism, a binary BEC with antiparallel dipolar components can be created. As shown in 
Ref. [21| and in this paper (see below), binary mixtures of the latter type may feature remarkable properties. 

The goal of this paper is to analyze the IMT in the binary dipolar BEC formed by the same atoms in different 
hyperfine states, in the presence of the linear coupling between the components. It may be expected that the 
nonlocal character of the dipole-dipole (DD) interactions may essentially affect the spatial separation between the 
components in the binary condensate. In accordance with this, the main issues we are addressing here arc effects of 
the DD interactions on the IMT, including a possibility to induce the IMT by the DD interactions without the linear 
intcrconversion between the constituents. 

The paper is structured as follows. In section II, equations describing the two-component dipolar BEC in the 
harmonic trap, including the linear coupling between the components, are introduced. This is a system of one- 
dimensional (ID) integro-differential equations, with the nonlocal DD interaction terms set down as in Ref. plj ]. 
Section III presents basic results obtained for the binary dipolar mixtures, with different combinations of orientations 
and magnitudes of the dipole moments in the two components. The effect of the DD interactions on the IMT is 
studied numerically, for a fixed total number of atoms, N\ + N2 = const. Dynamical properties of localized mixed 
and separated states are considered too. The paper is concluded by section IV. 



II. THE MODEL 



The quasi-lD dynamical equations for the dipolar BEC trapped in the harmonic-oscillator (HO) potential are 
derived starting from the rescaled Gross-Pitaevskii equation (CPE) with the short-range and DD interactions, 
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following Ref. 22]. In Eq. (TTJ), wave function ^/(r, t) is normalized to the total number of atoms, w p and uj z are 
normalized trapping frequencies in the transverse and axial (z) direction, p = \J x 2 + y 2 , 9 is the angle between the 
parallel dipoles and vector (^f—r'^j, and coefficients 7 and G measure, respectively, the strengths of the contact 
and DD interactions. The only assumption in the subsequent derivation is the standard single-mode approximation 
for the transversal part of the wave function, i.e., it is substituted by ^(z,p,t) — ip(z,t)4>o(p), where <fio(p) = 
(0FO -1 exp(-p 2 /2J 2 ), and / = 1/^u^is the characteristic harmonic-oscillator length which determines the transverse 
localization. The averaging in the transverse plane leads to an effective ID kernel of the DD interaction [2^ |. 



y DD (z) = ^[2^-^(i 



2r) exp (T)erfc(i/r)] 
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where t = {z/l) 2 , a = 1 and a = —0.5 for 9 — and 9 = tt/2, i.e., respectively, for the polarization of the dipoles 
along and perpendicular to axis z, and erfc is the complementary error function. Being interested in effects of the DD 
interactions on the IMT in binary BEC, we will focus on the basic situation, when the characteristic axial length scale 
is much larger than the transverse size, Z ^> I. In particular, the experimental setting with aspect ratio Z/l ~ 10 was 
realized in the chromium condensates [24I |25| . Under this condition, expression ^ yields results which are practically 
tantamount to those obtained with the use of the standard form of the DD interaction kernel, Vdd ~ 1/M 3 , truncated 
on the scale of z ~ I (the analysis performed in Ref. [23| for gap solitons in the quasi- ID condensate including the 
DD repulsion or attraction and a periodic potential, has demonstrated that the difference in the results generated by 
means of the different kernels is limited to < 3%). Therefore in the numerical calculations the DD interaction kernel 
was taken in the usual regularized (truncated) form. 

Proceeding to the consideration of the effectively ID two-component trapped BEC, we use a system of nonlinearly 
coupled GPEs, which include the DD terms, for the wave functions of the two components, t/ji(z,t) and tp2{z,t). If 
they represent different hyperfine states of the same atom, the interconversion, induced by a resonant electromagnetic 
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wave, is accounted for by linear mixing terms [26|, |27|. The scaled form of the respective GPE system is 
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where Vo is the strength of the axial trap, 711, 722 and 712, 721 are, respectively, the coefficients accounting for 
intra-species and inter-species contact interaction, k is the linear coupling, and A/i is a possible chemical-potential 
difference between the two components, which may be induced by an external dc magnetic field interacting with the 
atomic spins. Further, the DD interaction coefficients, Gu, G12, G21, G22, are defined as 



{Gii,Gi2,G 2 iG 2 2} = {vn\,{l- y)n\H2, i^i/x 2 , (1 - ^2} 
x (^ /4tt) (f - 3 cos 2 6), 
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where /io is the magnetic permeability of vacuum, /ii, [12 are the dipole moments in two BEC components, v = Ni/N 
denotes the ratio of the atom numbers in them, and N = Ni + N2 is the total number of atoms. The angle between 
the system's axis and the dipole orientations is 9, as before. Below, we focus on the two above-mentioned basic 
orientations: parallel (0 = 0) and perpendicular (8 — 7r/2), which correspond, respectively to the attractive and 
repulsive DD interactions. Two dynamical invariants of Eqs. © are the total norm, 
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and energy (Hamiltonian) 
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Stationary solutions are looked for as ipi^(z,t) = fa t 2(z) exp (— i/ii), where n is the chemical potential (in the 
presence of the linear coupling, the chemical potential of both components must be the same [l6j). and real functions 
obey the following integro-differential equations: 
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Properties of the binary BEC described by this model depend on a large number of parameters, including the 
number of particles in each component and the magnitudes and orientations of the dipoles. In the following, the 
results are presented for 711 = 722 = 10, 712 = 721 = 20, Afi = 0, Vq — 0.5, and equal numbers of atoms in both 
components, i.e., v — 0.5. As for the DD terms, we assume that the dipole moments of the two components are 
oriented in the same direction, with all the DD strengths being equal, Gn = G22 = G12 = G21 = T (case I), or with 
different strengths, Gn = T, G12 = G21 = LIT, and G22 = 1.21T (case II). As seen from Eqs. (j4|), cases I and II 
correspond, respectively, to equal magnitudes of the dipole moments, /ii = and unequal ones, (12 = The case 

of the BEC composed of two identical components polarized in opposite directions, Gn = G22 = —G12 = —G21 = T 
(case III) is considered too. The latter one corresponds to the above-mentioned setting composed of type-A molecules, 
in terms of the Hund nomenclature, under the action of the magnetic field, which determines the orientation of the 
magnetic moments, while the electric dipolar moments can be parallel or antiparallel to the direction of the applied 
magnetic field (2lj . 

In the absence of the DD interactions and linear coup ling , Eqs. <j3j) or ([7]) can be solved, under specific conditions, 
analytically by means the inverse scattering technique [26(, and numerically in the general case [28J. In particular, 
both the analytical and numerical solutions produce solitons. On the other hand, the localized miscible and immiscible 
states were found in the condensate with two components coupled by the linear interconversion. The IMT in the 
binary condensate can be controlled by the strength of the linear driving Hi . The critical value of the linear-coupling 
parameter for the IMT was estimated, by means of the variational approximation, to be n CI = (712 — 7n)/2, which has 
been confirmed by the numerical simulations. In this work, we use numerical methods to find fundamental solutions 
of Eq. ([3]) for the binary dipolar BEC, and explore effects of the DD interactions on the IMT in this system. 

III. IMMISCIBILITY-MISCIBILITY TRANSITIONS 

In this section, we present results corresponding to the fixed norm (rescaled number of atoms), P ~ 4, a free 
parameter being the linear coupling, n. As mentioned above, the parameters of contact interactions and trapping 
potential are the same as in Ref. 711 = 722 = 712/2 = 721/2 = 10 and Vq = 0.5, which correspond to the strong 
immiscibility, in the absence of the linear coupling and DD interactions. This may be realized so that the two species 
do not mix because their mutual repulsion is stronger than the intrinsic repulsion between atoms belonging to the 
same species. 

As concerns the parameters of the DD interactions, we will consider the following characteristic cases, I, II, and 
III (as mentioned above): Gn = G 12 = G 2 i = G 22 = T, or Gn = T, G 12 = G 2 i = LIT, G 22 = 1.2ir, or 
Gn = — G12 = — G21 = G22 = E Recall that T < and T > correspond to the attractive and repulsive DD 
interactions, respectively. Stationary states were obtained by simulations of Eqs. © in imaginary time. In the 
absence of the DD interactions, slightly perturbed stationary states maintain themselves (in real time) as breathing 
localized modes with the conserved norm 161. 



1. Equal dipole moments in both components 

In the absence of the linear coupling between the components, or with a very weak linear coupling, the addition of 
the arbitrary strong repulsive DD interaction to an initially immiscible state of the binary system cannot induce a 
transition to miscibility, see Fig. [TJa). However, the increase of the repulsive DD-interaction strength in the presence 
of a somewhat stronger linear c oup ling, which still corresponds to the immiscible phase in the case with only contact 
interaction, k < k ci (see Ref. [16[), can stimulate the IMT, as shown inQJb). In the presence of the repulsive DD 
interactions, the miscible state generated by the linear interconversion remains miscible, but the nonlocal interaction 
makes the amplitude of the trapped mixed state smaller, and its width larger, as seen in Fig. [He). These results may 
be summarized by stating that the repulsive DD interaction tends to mix the separated components, but it can do 
this only in the combination with the linear interconversion between them, while the interplay of the DD repulsion 
with the contact nonlinear interactions is less important, in terms of the qualitative description. 

The stability of the localized states was checked by direct simulations of their perturbed evolution in the real time. 
It was found that both the miscible and immiscible states are stable, although the former ones are more sensitive to 
asymmetric perturbations than their immiscible counterparts with the same norm. Generally, the perturbations give 
rise to small persistent intrinsic vibrations in the trapped states (not shown here). 

In the presence of the attractive DD interaction, the components that were separated at k = 0.1, in the absence of 
the DD interactions, remain separated. Naturally, with the increase of the DD attraction strength, the localized states 
become narrower and taller, see Fig. [5] However, localized states cannot be found for an arbitrarily large strength of 
the DD attraction (an explanation to this fact is given below). In other words, there is an upper limit of the strength, 
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FIG. 1: (Color online) The amplitude profiles of stationary states in the binary BEC for different values of the linear-coupling 
constant: (a) k — 0.1, (b) k = 2, and (c) k = 3. In each case, the profiles are displayed for the varying strength of the repulsive 
DD interaction, Gn = G22 = G21 = G\i = T. In this and other figures, parameters -y mn (m,n = 1,2), which characterize the 
contact repulsion, are 2711 = 2722 = 712 = 721 = 20. 




FIG. 2: (Color online) The amplitude profiles of stationary states of the binary BEC for different values of k: (a) k = 0.1, (b) 
K = 2, and (c) « = 3, in the presence of the attractive DD interaction with V = Gn = G22 = G21 ~ G12. The values of V < 
are indicated in the plots. 

beyond which the attractive DD forces cannot support localized modes. The latter conclusion is consistent with the 
real-time simulations, which confirm that the BEC components stay separated, unless the DD attraction is so strong 
that it destroys both components. An example of the evolution of the initially separated state in the presence of the 
strong attractive DD interaction, but yet below the destruction threshold, is shown in Fig. [3] for F = —10. Although 
the central part of the mode stays localized, the mode radiates a significant part of its norm. 

In the parameter region near the IMT induced by the linear coupling in the absence of the DD interactions [for 
example, at k = 2, see Fig. HJb)], the increase of the strength of the attractive DD interaction again makes the 
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FIG. 3: (Color online) The amplitude profiles versus time for a slightly perturbed immiscible state with n = 2, T = —10. Panels 
(a) and (b) display the evolution of the two components. The initial (dotted and dotted-dashed) profile and the one at t — 7.5 
(solid and dashed curves) are shown in (c). 
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FIG. 4: (Color online) The amplitude profiles versus time for a slightly perturbed mixed state with « = 2, T = —20. Panels (a) 
and (b) show the evolution of the two components. The initial very narrow profile (dotted black curves) and the one at t — 7.5 
(continuous blue and red curves) are plotted in (c). 




FIG. 5: (Color online) The blue lines with circles show the immiscibility-miscibility transition curves in the parameter space, 
(r, k), produced by the numerical calculations for cases I (co-polarized components, panel a) and III (counter-polarized 
components, panel b). The areas of miscible and immiscible states (above and below the transition lines with symbols) are 
labeled accordingly. Below the red thick line, localized states in the binary BEC with co-polarized dipolar components (case I) 
were not found. 



immiscible localized mode narrower and taller. Very strong DD attraction may result in a miscible state with a very 
high amplitude and very narrow shape. However, such states are unstable and disintegrate with time. An example of 
the instability is displayed in Fig. @] 

In the region of k slightly above K cr , i.e., in the states which are miscible in the absence of the DD interactions, 
the presence of the weak DD attraction may lead to a transition back to the immiscible state, which is unstable. In 
the rest of the region with k > K cr , the increase of the attractive DD interaction leads to taller and narrower miscible 
localized states, see Fig. [U(c). These newly formed miscible states are also unstable to small perturbations. Localized 
states cannot be formed if the DD attraction is too strong, similar to what was said above for smaller k. 

The nonexistence of localized states in binary BEC for arbitrary values of k in the presence of very strong attractive 
DD interactions is a consequence of the singularity in the ID condensate model, which manifests itself when the 
characteristic length of the DD interaction becomes comparable to the inter-particle distance [2l|. This is related 
to the singularity in the single-component BEC, in the case when the negative pressure caused by the inter-particle 
attraction is not compensated by the quantum pressure in the trapping potential (l9j . (2lj . 

Thus, the DD interactions can effect the IMT in the dipolar BEC mixture, with the components polarized in the 
same direction (and with equal magnitudes of the magnetic moments), which are coupled by the linear interconversion. 
This is a direct consequence of the nonlocality of the DD interactions. However, the DD interactions of either sign 
(repulsive or attractive), by themselves, cannot induce the IMT, as seen in Fig. [Sfa). On the other hand, the repulsive 
and attractive DD interactions' shift of the equilibrium in favor of the miscibility and immiscibility, respectively. 

A transition to the miscibility may be induced by the DD interactions (in the absence the linear interconversion) 
if the two dipolar components are polarized in the opposite directions, which was defined above as case III, with 
/ii = —fi2 in Eqs. ([4]). Recall that such a situation can be realized for molecules of the Hund's type A in the external 
magnetic field [21] . The possibility to induce the IMT by the DD interaction in this case is illustrated by Fig. H](a), 
which corresponds to the model without the linear coupling. 

The numerically found IMT curve for case III (the counter-polarized components) in the parameter space of (T, k) is 
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FIG. 6: (Color online) Amplitude profiles of stationary states for the two-component dipolar mixture with counter-polarized 
components: (a) in the absence of the linear coupling (k = 0) and (b) for k = 1. Parameters of the DD interaction are 
Gn = G22 = r, G21 = G12 = — r, with values of V indicated in the plot. The transition to the miscibility with the increase of 
r is evident. 
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FIG. 7: 2D plots of the evolution in case III (counter-polarized components) with k = 1: (a) and (b) the evolution of the two 
components initialized in the immiscible state for F = 2; (c) the evolution of a stable miscible state for T — 3.25, cf. Fig. 
[BJb). The immiscible states evolve as localized structures with complex intrinsic patterns, while the miscible state evolves as a 
regular breather. 

displayed in Fig. EJb). In this case, the addition of the linear interconversion shifts the transition to the miscibility to 
a lower strength of the DD interaction. Figure [7] displays examples of the evolution of immiscible and miscible states 
for k — 1 which are presented in Fig. [S^b) . The immiscible state evolves as a localized mode with a complex intrinsic 
structure, while the miscible state is only slightly perturbed, evolving into a regular breather. Note that, in case III, 
with the opposite polarizations of the dipoles in the two components, i.e., Gn = G22 = T < 0, G12 = G21 = T > 0, the 
IMT to a stable miscible state cannot be induced by the DD interactions in the absence of the linear interconversion, 
see Fig. [5] (b). In this case, only a narrow strongly pinned miscible state is formed, which is unstable. 

In fact, the configuration with the counter-polarized components considered here implies that, in both of them, 
the moments are directed along the spatial axis. The other case, with the counter-polarized components oriented 
perpendicular to the axis, is possible too, in which case the DD interaction is attractive between the components and 
repulsive inside each component. 

To summarize the analysis of cases I and III, we have found that the repulsive DD interactions affect the IMT, 
but only in the presence of the linear coupling in the case I- namely, the curve with circles in Fig. [SJa) shows the 
critical values of k which are lower than their counterparts found in the absence of the DD interactions. In this case, 
both the miscible and immiscible stationary modes remain stable under the action of the DD repulsion. The weak 
attractive DD interactions may convert a miscible state near the IMT threshold (k rs K cr ) into an unstable immiscible 
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one. At larger values of n > k ci , the increase of the attractive DD interaction causes the formation of tightly pinned 
miscible states, which are strongly unstable. There is a limiting value of the strength of the DD attraction, above 
which localized states cannot exist in the binary BEC. Finally, for the condensate with counter-polarized components 
(case III), the numerical results reveal the existence of the single IMT point, T, which moves to smaller values with 
the increase of k, see Fig. [5] (b). 

2. Unequal magnitudes of the dipole moments in the two components 

In the model describing binary BEC formed of atoms with equal dipole moments, which was considered above, it is 
obvious that all the states are spatially symmetric, in terms of the total density of both species (before and after the 
IMT, and for both the co- and counter-polarized components). On the other hand, in the binary system with different 
dipole moments of the two s pec ies an issue is not only the IMT itself, but also the equilibrium partition of atoms 
between the species [cf. Ref. [16[ , where an asymmetry between the species was introduced through the difference in 
the chemical potentials, Afi in Eqs. ([3]). 

Here we present results for slightly different magnitudes of the dipole moments, viz., Gn = r,G*22 = 1.21T, G 2 i = 
G12 = LIT, and several different values of T. In the presence of the repulsive DD interactions at small k, immiscible 
states are formed, in which the shapes of the two components may be spatially asymmetric. Depending on the 
initial profile used in the imaginary-time simulations, which were taken as two Gaussians with coinciding or separated 
centers, two distinct types of stationary immiscible states could be found, respectively, in the binary system: those 
which feature the mirror symmetry between </>i(z) and foiz), and asymmetric ones, see Fig. [8] Computing energies of 
the observed states, we have found that the energy is smaller (by a few per cent) in the asymmetric state. Nevertheless, 
it does not represent the ground state of the system, because in direct simulations this stationary state develops into 
a breather, as shown in Fig. [SJ Note that the component with the larger dipole moment features a smaller amplitude. 
Opposite to the case of equal dipole moments, where the growth of the repulsive DD interaction could lead to the 
IMT at k ^ 0, here the growth of the DD repulsion increases the separation between the components. A noteworthy 
difference of the immiscible states in the present case from their counterparts in the model with equal dipole moments 
(cf. Figs. [U [21 O for the case of the co-polarized components, and Fig. [5] for the counter-polarization) is that 
the unequal moments give rise to three-lobe shapes. They feature two regions occupied by the component with the 
larger dipole moment, which are separated by a single region in which the component with the smaller moment is 
concentrated. 

The character of the mixed states generated due to the linear coupling is not changed in the presence of the 
repulsive DD (of course, the shapes of the two components become different). Real-time simulations indicate that 
the immiscible and miscible localized states in the binary system with different dipole moments of the components 
evolve into breathing structures, which tend to keep the original immiscible or miscible structure of the pattern. This 
is the case for both symmetric and asymmetric profiles of the two components, which is understandable due to very 
close energies of the immiscible and miscible modes (the energy difference between them is < 1%). An example of 
the breather developing from immiscible modes is displayed in Figs. |FJ(c,d). 

The attractive DD interaction pushes all binary states in the system with unequal dipole moments to enhanced 
immiscibility, see Fig. [9] (which is consistent with effects of the DD attraction outlined above in the symmetric system) . 
As well as in the case of the DD repulsion, the symmetry of the stationary states obtained from the imaginary-time 
simulations depends on the initialization, cf. Figs. 151 and HUT aLfb). The separated asymmetric states are getting 
more and more asymmetric and both symmetric and asymmetric ones are becoming more and more narrow with the 
growth of |r|. The maximum energy difference between the symmetric and asymmetric states for fixed T is of the 
order of a few per cents. The stationary mixed state which exists without the DD interaction for certain values of k 
is transformed by the DD attraction into immiscible states whose symmetry is determined by the initial conditions, 
see Fig. HUT b). Note that, in the case of the DD attraction, the stationary localized states can only be found if the 
DD interaction strength is not too large. Real-time simulations demonstrate that all localized states which exist in 
the presence of the attractive DD interaction are unstable, see Figs. ITU1 (c)-(f). 

The results presented above, and additional ones not displayed here, that were obtained for different combinations 
of the DD interaction parameters, show a great diversity of possible localized states in binary BEC with components 
formed of atom with different moments. Generally, their symmetry properties, shape, and width depend on the initial 
conditions. 
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FIG. 8: (Color online) Profiles of the trapped modes in the binary BEC with k = 0.1, produced by the imaginary-time 
simulations which were initialized by separated (a) or overlapping (b) Gaussians. In the absence of the DD interactions, the 
corresponding state is immiscible, independently of the initialization. In the presence the DD interaction, the modes evolve 
into breathers, as plotted for both components in (c) and (d) for T = 40 [the respective initial profile is presented by the dashed 
line in (b)]. Parameters of the DD interactions are Gn = T = 40, G21 = G12 = LIT, and G22 = 1.21D 




FIG. 9: (Color online) Profiles of asymmetric stationary states in the binary BEC for different values of the linear-coupling 
constant: (a) k = 0.1, (b) k = 2, and (c) k — 3. The strengths of the attractive DD interaction are Gn = F, G22 = 1.21F, and 
G21 = G12 = l.ir. The respective values of T are indicated in the plots. 
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FIG. 10: (Color online) Profiles of the localized modes in binary BEC with k = 3, obtained by means of the imaginary-time 
simulations with initial conditions in the form of separated (a) and overlapping (b) Gaussians: T = (solid lines), P = 40 
(dashed lines), and V — — 10 (dotted lines). In the absence of the DD interaction, the corresponding state is miscible, irrespective 
of the initialization. The states are unstable for either sign of the DD interaction, as shown in (c), (d) for r = 40 [the respective 
initial profile is shown by the dashed line in (b)] and in (e), (f) for T — — 10 [the initial profile is presented by the dotted line 
in (b)]. In both cases, the twin panels display the evolution of the two components. Parameters of the DD interactions are 
Gn = T, G21 = G12 = l.ir, and G22 = i.2ir. 
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IV. CONCLUSION 

In this work, wc aimed to study effects of the DD interactions on IMT in binary BEC trapped in the HO potential. 
The interest to this problem originates from the possibility to study the interplay of the nonlocal DD attraction or 
repulsion with the global spatial structure induced by the IMT. This binary system is a complex one, as it depends on 
many parameters, such as the orientation of the dipoles and the strength of the interaction between them, numbers of 
atoms in each component, and the rate of the linear interconversion between the components induced by the external 
spin-flopping electromagnetic wave, in case the two components represent two different spin states of the same atom. 

We have reported the results for three dipolar binary BEC systems: with co-polarized or counter-polarized dipolar 
components and equal magnitudes of the dipole moments in the two components, and for the system with unequal 
moments. Both repulsive and attractive DD interactions were considered. In all the cases, the IMT crucially depends 
on the presence of the linear coupling between the two species, while the location of the miscibility threshold and 
the stability of the resultant immiscible and miscible states may be significantly affected by the DD interactions. In 
general, the repulsive and attractive DD interactions shift the equilibrium in favor of the miscibility and immiscibility, 
respectively (i.e., the long-range repulsion tends to suppress the spatial structure induced by the immiscibility, while 
the nonlocal attraction helps to enhance it, although with the trend to make it unstable against the spontaneous 
transformation into a breather). In those cases when regular breathers develop from unstable immiscible or miscible 
stationary modes, they tend to keep the original immiscible/miscible structure of the pattern. Only in the system 
formed by the components with equal counter-polarized dipole moments, the IMT may be induced by the DD interac- 
tions in the absence of the linear interconversion between the constituents in the binary BEC. On the other hand, in 
the case of the attractive DD interactions, stable localized states do not exist for too strong attraction. In particular, 
slightly perturbed miscible narrow states, formed in the binary BEC with a very strong DD attraction, are highly 
unstable and quickly disintegrate. 

The patterns predicted in this analysis may be realized by means of the available experimental techniques. In 
particular, the chromium condensate may be created in a quasi- ID trap with the aspect ratio ~ 10, the axial length 
about 10 /im, and the total number of atoms ~ 25,000 (2J,[25j]. This atom number and length should be sufficient 
to allow the observation of the spatial structures expected in the immiscible patterns considered above. The actual 
strength of the DD interactions between the chromium atoms is equivalent to the effective scattering length ~ 0.75 
nm, which can be made comparable to the strength of the contact interactions if they are properly attenuated by 
means of the Feshbach- resonance technique [24| . 

It may be interesting to extend the analysis reported in this work to a 2D setting. It should also be relevant to 
extend the model by taking into regard additional physical factors which are known to occur in dipole condensates, 
such as atomic loss due to inelastic collisions [25| . 
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